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TECHNICAL NOTE 2702 


AN APPROXIMATE METHOD OF DETERMINING THE SUBSONIC FLOW IN AN 
ARBITRARY STREAM FILAMENT CF REVOLUTION CUT BY 
ARBITRARY TURBQMACHIHE BLADES 
By Chung -Huh Wu, Curtis A. Brown, and Vasily D. Prian 


SUMMARY 

A method is presented to obtain a relatively quick approximate 
determination of the detailed subsonic flow of a nonviscous fluid past 
arbitrary turbomachine blades. The governing equations are formulated 
as thou^ the flow were constrained to stream filaments of revolution 
(the thicknesses of which vary as they pass throu^ the machine) . 
Attention is first fixed on a certain streamline lying between the 
blades and in such a filament. The shape of this streamline and the 
specific mass flow (product of density and velocity) along it could be 
estimated, for example, from the shape and the width of the channel 
formed by tiro adjacent blades and from the inlet conditions. From 
these starting values, the governing' equations yield separate values 
of density and velocity components and their derivatives in the circum- 
ferential direction. These values are Incorporated in a Taylor’s series 
to provide the analytical continixatlon of the flow quantities from the 
chosen streamline to the blade surfaces. Considerations of integrated 
mass flow then show how the assumed shape of and specific mass flow 
along the chosen streamline must be adjusted to provide a starting point 
in a second cycle of computation. 

The method is illustrated with examples of compressible flow in a 
t\irbine cascade and in a centrifugal conipressor. For these high-solidity 
blades, three terms of the Taylor's series are found to -give sufficient 
accuracy. Sufficient convergence is obtained in the turbine cascade 
after two cycles of ccairputatlon (started with an available incompressible 
mean streamline) and in the centrifugal compressor after four cycles of 
computation (without any aid of available information) . The detailed 
flow variation obtained compeLres very well with an available numerical 
solution and experimental data for the turbine cascade and with detailed 
experimental measurements for the centrifugal compressor. Because each 
cycle of computation for compressible flow takes only 16 hours, suc- 
cessive Improvement of the solution is practical. 
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lUTRODUCTaiOIT 

A 'basic aerodynamic problem of gas-turbine engines is the flow of 
air or gas through the compressor or turbine bladings. "When the dis- 
tance between hub and casing is relatively short, the air or gas may be 
assumed, for an approximate solution, to flow along surfaces of revolu- 
tion which are usually noncylindrical (figs. 1 and 2). The equations 
governing the flow of a nonviscous cong)ressible fluid along such sur- 
faces have been formulated recently by the use of a set of orthogonal 
coordinates . I and cp (meridional and angular, respectively) on the 
surface (figs. 1 and 2), and general methods of solution for both sub- 
sonic and supersonic floif are given (references 1 and 2) . In the case 
of supersonic flow, solution by the use of the method of characteristics 
is satisfactory. But in the case of subsonic flow, the general numerical 
method suggested takes a long time to accomplish if a high-speed large- 
scale digital computing machine is unavailable. In such a case, a much 
quicker approximate method of solution which gives sufficient accuracy 
is desirable. 

■5-Jhen the distance between hub and casing of the cconpressor or tur- 
bine is relatively large, the two-dimensional solution of the flow along 
the surfaces of revolution cut by the blade sections becomes inadequate. 

A general theory of three-dimensional flow is given in reference 3, in 
which the coniplete three-dimensional flow in a tirrbomachine is obtained 
by a suitable ccmbination of t^fo-dimensional flows along two kinds of 
stream surface extending, respectively, from hub to casing and from 
blade to blade (fig. 3). In the first stage of such calculations, it 
is desirable to have some general approximate knowledge of the detailed 
flow variation from blade to blade obtainable by relatively quick 
approximate methods by assuming that the gas flows along surfaces of 
revolution. 

For these purposes, a quick approximate method of solution of the 
detailed subsonic flow of a nonviscous fluid past arbitrary turbomachlne 
blade sections along an arbitrary surface of revolution has been 
developed at the HACA Lewis laboratory. The method is essentially an 
extension of the general- two-dimensional blade design method reported 
in reference 4 and is particularly useful for high- solidity blades, such 
as those encountered in axieQ.-flo\f turbines and radial- or mixed-ficrt'r 
compressors and tirrblnes. The method is presented and is then illus- 
trated by two examples, namely, con 5 >ressible flow past a turbine cas- 
cade and compressible flow along a curved stream filament of revolution 
(at the mean blade height) in a centrifugal ccmipressor. The velocity 
and pressure distributions along the blade surfaces as well as across 
the channel are conqiared with available numerical solutions or experi- 
mental measurement & or both. 
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SYMBOLS 

The following symbols are used in this report: 
h static enthalpy 

I h + I W2 - I aj2r2 

L blade length projected on turhomachine axis 

Z,cp orthogonal coordinates on mean surface of revolution 

M mass flow 

N number of blades 

P pitch or spacing 

p static pressure of gas 

p^ total, pressure of gas 

r radial distance from axis of turhomachine 

T normal channel distance from huh to casing in meridional plane 
t hlade thickness in circumferential direction 

U hlade speed, cur 

V absolute velocity of gas 

W velocity of gas relative to hlade 

y distance equal to rcp for flow on cylindrical surface 

z distance along axis of turhcaaachine 

B slope of streamline on stream surface, tan"-^ r/- or tan"-*- ^ 

^z 

y ratio of specific heats 

A correction in specific mass flow 

p density of gas 

a slope of meridional curve, tan"-*- — 
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T normal thickness of stream, filament of revolution 

ijr stream function 

£0 angu3-ar velocity of "blade 

Suh scripts: 

c on chosen streamline 

1 at inlet 

I,cp meridional and circumferential components, respectively 
m on mean streamline 
p pressure surface of hlade 
s' suction surface of hlade 
t total or stagnation state 
y y- component 

2 z-component 
Superscript; 

* dimensionless' value obtained hy dividing hy inlet value 


METHOD 

Equations Governing Flow on Arbitrary Surface 
of Revolution 

■When the gas is assumed to be nonviscous and to flow steadily along 
surfaces of revolution in turbomachines, a relatively simple description 
foimulated in reference 1 qf the governing equations for the gas flow 
is applicable. For the gas flowing in an arbitrary stream filament of 
revolution defined by the meridional coordinate. I and the angular 
coordinate of the mean surface of revolution and by the varying 
normal thickness T of the stream filament (fig. l), the continuity 
eq-uation takes the following form: 

^(TpW^r). 5(TPWjp) 

3cp 


= 0 


( 1 ) 
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The circumferential campcaient of the equation of motion for irrotational 
absolute flo^r is 


. hVcp 1 hMli 

bl r. ^ 



sin 0 = 0 


( 2 ) 


The density of the gas is related to its value at the inlet and the 
local gas velocity hy the use of. the isentropic relation between enthalpy 
and density of perfect gases 


)* = - 2 - = 


I + i oi^r^ 





(3) 


in which I is a constant for steady adiabatic flow. 


General Description of Approx im ate Method of Solution 

In order to obtain a relatively quick approximate solution of the 
flow equations (l) to (3) on a given arbitrary stream filament of 
revolution cut by axbitrary turbomachine blades, attention is first 
fixed on a particular streamline about midway between two neighboring 
blades (fig. 2(a)), as in the design method of reference 4. The shape 
of this streamline and the specific mass flovr along it are esti- 

mated from the given blade shape and the inlet and exit conditions, and 
separate values of velocity ccmrponents and density on this chosen 
streamline are then easily computed. After this, successive clrcum- 
ferentleil derivatives of velocities are obtained from the values 
obtained on the chosen streamline, and the fluid properties in the 
angular direction are obtained by eniploying these derivatives in a 
Taylor’s series. To this point, the calculation is exactly the same 
as that in the blade design method given in reference 4. 

After completion of this phase of the work, the mass flow between 
the chosen streamline and the two given blade surfaces is computed, 
rather than, as in reference 4, the blade coordinates being computed for 
a ^ven inlet mass flow. By comparing the relative percentages and the 
total mass flow with the value given at the inlet, the phape of the 
chosen streamline and the specific mass flow along it are corrected. In 
general, this process is repeated tintll convergence is obtained, (in 
the subsequent examples, after seme experience on the part of the com- 
puter, the later cycles for both examples are ccmpleted within 16 hr 
each.) The acciaracy of the solution of the problem then depends most 
,on the number of terms used in the Taylor's series. For a check of the 
accuracy, the data obtained in the solution lead quickly to a calculation 
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of the residuals of a second-order partial differential equation in the 
stream function, which is obtained hy combining equations (l) and (2) 

( see reference 4) . If the residuals are found to be large, the solu- 
tion can be readily improved by the relaxation technique and a more 
accurate solution of the given blade on the given surface is established. 


Choice of Starting Streamline 

The choice of the staxtlng streamline used in the calculation 
depends on individual cases. In axial-flow-type machines, the stream- 
line circumferentially dividing the mass flow in the channel into ti^ro 
equal parts (fig. 2(a)) is best to use and is found to resemble some- 
what the blade mean line in reference 1. This streamline is designated 
the mean streamline, and the flow on it may be taken as the mean in the 
channel. In the centrifugal con 5 )ressor problem to be presented later, 
the mean streamline is found to be too close to the suction surface of 
one blade. (The shape is also quite different from the blade mean line.) 
The extension of the flow from the mean streamline to the pressure sur- 
face of the outer blade therefore requires more terms in the series. 

As a compromise, a streamline having 65 percent of the mass flow betvreen 
it and the suction surface of one blade is chosen. 


Calculation of Flow along Chosen Streamline 

If no better Information is available, the shape of the chosen 
streamline is taken as that of the blade mean line and the specific mass 
flow along it is estimated as the average value between the blades 
according to the following relation: 

■ = JL (4) 

(PWz)iTiri P-t ^ ^ 


or 


(P^l^c - IWiJI 


'^'i^i P 
Tr P-t 


(4a) 


(See the appendix for the scheme of nondimepsionalization. ) From this 
value of "the slope of the chosen streamline (tan p^,), and the 

given inlet condition, separate values of W^, and p along the 

chosen streaiallne are obtained as follows: Equation (3) is rewritten 


o 


2384 



MCA TN 2702 


7 


2 = 



where 


(5) 


0 = (p*W2 



( 6 ) 


and 


2= p! 


I + i CD^r^ 


1 

r-1 


(7) 


Tables listing 2 in small intervals of $ are given in refer- 
ence 3 for T = 1.4 and 4/3. (There is a misprint of a minus sign in 
the exponent of equation (5) at the top of the tables given in refer- 
ence 3. ) With the given values of p*W* sec p along the chosen stream- 

line^ $ is calculated according to equation (6), Z is then read from 
the table, and p* is obtained according to equation (7). After p* 

is determined, W? and W^ are easily obtained by use of the fol- 
’ 1 , 0 . cp^ c 

lovring formulas: 


W* = 

”z,c 


(pX). 


( 8 ) 


W 


= W? tan p 


9,c " 1,0 


(9) 


T-Jhenever there is available information concerning the flow ^last 
similar blades, this information should be used to obtain a better esti- 
mate of the shape of and the flow on the chosen streamline in the first 
calculation. 


Angular Variation of Fluid State 

The first- and second-order angular derivatives of W^, W^, and p 

at the chPsen streamline can be obtained from the basic equations (l), 
(2), and (3) (see reference 3 for derivation) as 
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Sep 




^cp •fcan 3 

d2 + m ■" '"'P 


+ (%) + 2asr) sin a 


coB^ p (lO) 


Sw, 


W 


5^ 


dW, 


”Q) 

r + (Wq) + 2mr) sin a 


■tan p - — 


2 _ d(TrpW^) 


Tp dZ 




> cos^ P (ll) 


1 Sp 
P 5^ 


/ S¥x h\kp 

rsr «i 5?- + '''P5v 


(r-1) (I + I - I T^?2) 


(12) 


_ J 2 Sp 1 d 

Sep^ I Tp2 ^ **"P ^ 


SW, 


tan P 


( A . 

dl dcp + 


sin CT 


dqp 


7 cos^ p 


(13) 


d % a^w 


(14) 


1 ^ _ 2^ ( Sp"\^ 

<>3,^ p2 'v^, 


(r-1) (1 + § “^r2 - i w2) 


S%7 /SW7\2 /S¥co\2 


Sqp' 


(15) 


Third- and hi^er-order derivatives, if reejuired, can also "be 
ohtained. The variation of any fluid property q across the c h a nn el 
at a constant value of I is then computed hy a Taylor's series in 

(cp - qPc) 


l(cp) = ^ + 


(cp - qp^) 


(si): 


(cp -cpj' 


21 


S2g\ (cp-epe)^ 
S<p2;c 31 



( 16 ) 
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After the variation of fluid state ‘In the angular direction is 
obtained^ the mass flow across any ccmstant-Z line from the chosen 
streamline to the suction and pressTire surfaces of the two adjacent 
blades is ccmputed by the following formulas (fig. 4): 

r^c 

Ms=TrJcpg dqp _ (l7) 

Mp = “Tr Jcp^ .pWl dcp (18) 


Correction of Shape of and Specific Mass Flow along 
Chosen Streamline 

The correctness of the solution within the accuracy of the finite 
terns retained in the series is indicated by the following three items; 

(1) the closeness of the smi of Mg and to the given inlet value, 

(2) the constancy of the relative magnitudes of Mg and Mp along the 

streamline, and (3) the closeness of the ratio of Wcp to con 5 >uted 

at the blade surface and the slope of the given blade siirface. 

A simple method of correction- is suggested. In figure 4(a) are 
shown a number of stations used in obtaining the mass flcrtf across a 
constant-Z line, which are equally spaced in each of the two segments. 

The variations of pW^ and mass flow along cp obtained in the solution 
are shown by the solid lines in figures 4(b) and 4(c), respectively. 
VJhether (Mg + Mp) is greater or smaller than the inlet mass flow indi- 
cates whether the p¥^ value is, in general, higher or loiter, respec- 
tively, than the correct veilue. With the variation of pW-j id.th respect 

to cp assvuned correct, a constant (with respect to cp) a is added to - 
the PW^ values obtained in the solution to fit the mass flowj thus 




A = 


Mi - Tr Jcp PWj dcp 

5 

Tr (q?p - cPg') 


(19) 


A corresponding increase in mass flow equal to Tr(cp-<¥>g) multiplied by 

A is then added to the curve in figure 4(c), resulting in the dashed 
line. A new position of the chosen streamline cp^ is obtained by read- 
ing from this new curve at the chosen percentage of mass flow. 
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This process of correction is to "be applied at each station along 
the mean streamline. Fiorther adjustment can he made for the right flow 
direction at the hlade surfaces. Q5bie chosen streamline may also he 
corrected a short distance outside the channel hy using the finite- 
difference equations given in reference 1, or end-point differentiation 
formulas may he used at the first and last stations inside the chemnel. 
If accurate determination of the velocity distrlhution near the leading 
edge is desired, either shorter spacings can he chosen there and dif- 
ferentiation coefficients for unequal spacings (reference 5) can he 
used, or a portion of the flow region near the leading edge can he 
in 5 )roved hy the finite-difference methods given in reference 1. 


AEPLIGAIIONS 

The method descrihed is applied to obtain an approximate solution 
of con 5 )ressihle flow past a turhine cascade of uniform T and along a 
curved stream filament of revolution -vriLth variable t in a centrifugal 
compressor. 


Turhine Cascade 

This is a relatively simple case of the general method. For two- 
dimensional flow past a cascade of blades, T is uniform throughout and 
r/ti is unity; also, r is a constant and sin a = 0. For convenience, 
the coordinates r and q) are combined to become y (coordinate I 
becomes coordinate z) . 

The cascade chosen is that investigated in reference 1, for which 
data shoTfing detailed theoretical flow variation at inlet Mach numbers 
of 0.42 and 0 and the e3g)erimentally determined pressure distribution 
on the blade at an inlet Mach number of 0.42 are available. 

The shape of the mean streamline and the variation of W* obtained 
in the lncon 5 )resslble solution (see figs. 8 and 10 of reference l) are 
■used as the starting mean streamline and the variation of specific mass 
flow p*W* along the mean streamline, respectively, for the compressible 
solution at an inlet Mach number of 0.42. 

After the flow on the mean streamline is determined, the extension 
in the y-directlon is accomplished by the use of the first three terms 
of the Taylor’s series. The mass flow computed for the given blade 
shape is found to b.e somewhat high in most places and the division of 
mass flow is such that Mg is slightly lower than Mp. A correction is 

then made according to the method described. The results of the second 
solution are entirely satisfactory and no fxarther correction is made. 

(The computation follows the setup given in reference 4.) 
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The starting and corrected mean streamline coordinates are sh.own 
in figure 5(a) . The starting and corrected values of and the 

ratio P/(P-t) are shown in figure 5(h). The variation across the 
channel of the magnitude of the dimensionless velocity components as 
well as of the resultant velocity is compared with the relaxation solu- 
tion of reference 4 in figure 6. The velocity distrlhutlon around the 
blade is cor 5 )ared with the relaxation solution and experimental value 
'of reference 1 in figure 7. 

t 

In general, the compeirison is satisfactory. The good correlation 
is especially apparent for most of the* channel along the pressure sur- 
face of the blade. However, the velocity gradient of the entrance and 
exit sections for the suction siirface of the blade is somewhat dif- 
ferent; the experimental curve has a steeper gradient than either 
theoretical curve. 


Plai'f along Ctirved Stream Filament of Revolution at Mean Blade 
, Height in Centrifugal Inpeller 

The method is applied to a more general case of the flow along a 
curved surface of revolution in a centrifugal impeller for which the 
e 25 >erimentally measured blade-to-blade flow variation at the mean blade 
height is available (references 6 and 7). In reference 7, the theo- 
retical velocity distribution around the blade surface is obtained 
essentially by assroning that the mean streamline has the same shape as 
the blade mean lin e and that a linear variation of air velocity exists 
across the passage. Negative theoretical velocity, which is not observed 
experimentally, is obtained at the pressure surface . With the cor- 
rection procedure for the shape of and the specific mass flow along the 
chosen- streamline and the provision of a Taylor's series to consider 
more than a linear variation of the flow condition across the passage, 
the present method is applied to this inpeller for a more complete com- 
parison with the experimental data. 

Simplifying assunptions . - So that this theoretical computation 
can be made and the results compared with the experimental data, the 
foUovring assumptions must be emphasized; 

(l) In order to compare the theoretical results iriLth the available 
experimental data measiored along the surface of revolution at the geo- 
metrical mean blade hei^t (jis* 8), it is assumed here that the air 
peirtlcles lying on a circle of the mean radius 8.15 at the inlet sub- 
sequently form a stream surface of revolution which coincides with the 
geometrical mid-blade-helght surface instrumented in the eiforementioned 
experimental study. In general, the stream surface so formed by the 
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air particles lying originally on a circle at the inlet deviates from a 
surface of revolution. The magnitude of the deviation^ however, is yet 
to he determined. A surface of revolution may he assumed for simpli- 
city; if desirable, however, the position of the surface may he more 
accurately obtained from a huh -to-casing calculation along an S 2 
surface similar to that shown in figure 3 by the method given in refer- 
ences 3, 8, 9, or 10. This is not done herein because experimental 
data are available only along the mean-blade -height surface of revolu- 
tion. 


( 2 ) In the present solution, it is assumed that the ratio of the 

normal thiclmess of the stream filament T to its inlet value is the 
same as the ratio of the normal distance between hub and casing T at 
the same place to that at the inlet; that is, is equal to T/Tj^. 

A more correct value could again be obtained from the aforementioned 
ceilculation. From the results obtained in an incompressible solution 
given in reference 9 for a similar centrifugal impeller, it is found 
that f/Tj^ is larger than T/T^ in the first half of the flow path 
and sma ll er in the second half (fig. 9) . G?his approximation will there- 
fore give lower and higher velocities, respectively, for these t\To 
portions. 

( 3 ) Viscosity effect is neglected in the theoretical calculation. 

The boundary layer in the actual flow reduces the effective channel area 
and thereby increases flow velocity, and the viscous losses increase the 
entropy of the air and reduce the rise in static pressure. These effects 
should be small at the inlet and Increase gradually tc^'rard the exit. 

Computation procedure . - Without any help of available information 
on the relation between the shapes of the mean streamline and the mean 
blade line for this type of flow (on a general surface of revolution), 
the starting calculation is made with an assumed mean streamline the 
same as the blade mean line and located midway in the channel. Its posi- 
tion downstream of the blade is determined on the basis of a slip factor 
of 0.91 and the condition that the angular momentum of the air remains 
constant eO-ong the streamline. The starting value of is based 

on equation (4a) iTith an inclusion of a multiplication factor of 1.23 
to account for the reduction of channel cross-sectional area caused by 
the taper in the blade section ( see fig. 3 of reference 7 ) . Eleven sta- 
tions 1.1 inches apart are used Inside the impeller (figs. 2(b) and 8). 
The determinations of the flow along the mean streamline and of the first 
and second derivatives of flow velocities in the angular direction are 
given in table I. The determination of the mass flow between the mean 
streamline and the two blade s\irfaces at one station is show in table II. 
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The resiilts oltalned In this calculation show tliat the mass flows 
across the channel at most . of the stations are lower than the given 
inlet value and that the division of mass flow is such that the mean 
streamline should he shifted considerably toward the suction surface of 
the blade, especially along the middle portion of the streamline. So 
that the location of the chosen streamline mil not be too close to the 
suction side, a streamline which has a mass flow between it and the 
suction side of 65 percent of the total mass flow in the channel is uaed 
for the remaining calculations (fig. lO) . The variations of (**^*^*P*^'^)c 

and taken from the corrected values are shovm in figures ll(a) 

and U(b), respectively, by the square symbols. 

With these starting vadues for the 65 percent streamline, a new 
cycle of calculation is made according to the setup shown in tables I 
and II. Correction in the shape of and bhe specific mass flow along 
the 65 percent streamline is further made in the same manner as pre- 
viously. 

The same process is repeated t\ro more times. To expedite conver- 
gence, the corrected values obtained are not directly used as the start- 
ing value for the next cycle. Instead, an Interpolated value of all 
previous calciilations is used. This interpolation is done by plotting 
the starting and corrected values of each cycle of calculation along the 
abclssa and ordinate, respectively. The intersection of the curve join- 
ing these points and a 45° line passing throu^ the origin is taken as 
the more correct value to be used as the starting value for the next 
cycle of calculation. The third cycle of calculation thus obtained for 
the 65 percent streamline shows sufficient convergence. The successive 
shapes of T*r*p*W* and P*^J^ for the 65 percent streamline are com- 
pared in figures 11(a) and U(b) . The accuracy of the last solution 
is also indicated by the values of residuals (table Hi) of the differ- 
ential equation in the stream function computed according to the method 
given in reference 4. The estimated error in P*W^ divided by the 

approximate value is given in table IV. (These ratios indicate the 
order of magnitude of error in velocity. ) 

Comparison of theoretical and experimental flow variation . - The 
velocity distributions along the 10 percent and 90 percent geometrical 
lines taken from the experimental data reported in references 6 and 7 
are compared id-th the theoretical value computed for the same position 
in figure 12. It is seen from this figure that the theoretical velocity 
gives a qualitative pictiire of the experimental variation, and that the 
average theoretical velocity is, in general, relatively low, especially 
in the first half of the flow passage. If the effect of the assumed 
surface of revolution for the actual twisted surface is disregarded, 
most of this difference p 2 robably can be attributed to the relatively 
large error in the assumed T variation for the first half , as 
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indicated by figure 9. Toward the trailing edge of the blade, the 
reduction of flow area due to bo undar y layer along the blade surface 
should also increase the velocity in the actual case. It is of interest 
to note here that no negative velocity occurred at the pressure siirface 
of the blade in the theoretical solution. 

The theoretical and ejcperimental pressiires along the blade surface 
are compared in figure 13. In general, the two agree much better than 
in the case of velocity distribution, and the difference in pressure 
on the t^fo sides of the blade for the experimental data is larger than 
that for the theoretical value. The differences among the theoretical 
and experimental values should be explained through the difference in 
velocity and the smaller increase of total pressure due to losses in the 
actual, flow. 

The theoretical and experimental variations of relative velocity 
and static pressure across the passage are compared in figures 14 and 15. 
In figure 14, the theoretical velocity is lower than the experimental 
value in most places. From pressure surface to suction surface, the 
theoretical velocity also has a higher gradient than that of an average 
line drawn through the experimental data. Besides the assumptions 
made in the theoretical calculations, there is also the experimental 
error in the measurements of static and total pressure; this differ- 
ence in pressures is used in the computation of the experimental 
velocity. VThen the velocity is small, the normal error in the pres- 
sure measurements can cause a very large error in the velocity so 
computed (see reference 7) . 

On the other hand, the theoretical and experimental static pressures 
agree well throvighout the passage (fig. 15) . This good agreement indi- 
cates that the present approximate method of solution should be useful in 
predicting, comparing, or improving the performance of, at least, slmlleir 
type machines . 


CONCIUDING REMARKS 

A method is presented to obtain a relatively quick approximate 
solution of the detailed subsonic flow of a nonviscous fluid past arbi- 
trary turbomachine blades along a given arbitrary surface of revolution. 
The method is illustrated with examples of compressible flow in a tur- 
bine cascade and in a centrifugal compressor. Three terms in the series 
to account for the circumferential flo\7 variation from the chosen stream- 
line about laidway in the channel to the blade surface are found to be 
sirfficient. For a centrifugal compressor, the mean streamline shape is 
found quite different from that of the blade mean line; consequently, 
correction for the mean streamline shape is necessary. Sufficient con- 
vergence is obtained after two cycles of conputatlon for the turbine 
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cascade and four cycles for tlie centrifugal compressor. (Each, cycle of 
congnitatlon can he done In 16 hr. ) The detailed flow variations ' 
obtained compared very well with an available numerical solution and 
experimental measurements. 


Lewis- Flight Propulsion Laboratory 

national Advisory Committee for, Aeronautics 
Cleveland, Ohio, Februeiry 12, 1952 
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APjEEWDIX - rf0M)IMElNSI0MLI2ATI0N OF EQUATIONS 
In the exaiEples, all quantities are made dimensionless as follows; 


♦ r 

r = 


. T* = ^ - Z* = i- • W? = ^ • 

^i ’ '^1 ’ ^ ^l,± * ^ ^1,1 ' 




W* = 




• u* = - • p* = -fi- . h* = 

'h,± Pi w: 


^ ; I* = ; M* = 

fT ^ '^±^±Pih,± 

V ^ X 6 ^ X 


In this dimensionless form, equations (3), (4a), (6), (7), (lO) 
to (14), and (16) to (18) appear as follows: 


P* = 


'l* + i i ¥*V-1 




(3') 


jc “ P-t 


(4a* ) 


$ = 


tj*2 

2 '"'Z.i 


(pX sec p); 


'I* + I uf 


( 6 >) 


S = P' 


I* + I uf r*^^ 

~~K 


(7-) 


5^ 


tan R <i(T^r*P*VJ*) f . .A 

^ ^ — h- + wl+2ufr* 

dZ T*p* dZ* V ^ ^ / 


sin a 


cos^ P (lO') 
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(^5 + 2U*r*j sin 


tan p - 


1 d(T*r*p*wJ‘) 


T*P* dZ* 


y COS^ p 
( 11 * ) 


1 So* 


(r-1) (I* + i U|^ r*^ - i ¥*2) 




( 12 ') 


2 ap* *i(T*r*p*wf) 




^♦p*2 5^ dT*" 




T*p* 





,* d H 




tan p \r* sin a^^l^cos^ P 








+ 61n a 


S¥( 


^ cp 

sr 


(13*) 


(14') 


«? = K)c -'po) + lr(5^) 


,)^ + . . . (16*) 


M* = T*r* Jcp^‘= p*W* dtp 

(17') 


(18') 


In the second example, the velocity and static pressure are given 
in terms of, respectively, the velocity of sound at the inlet total con- 
dition and the inlet total pressure 
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p ^ _EL. 

®t,i ^t,i 

The dimensionless form of the partial differential eqnation in the 
stream function used to con 5 )ute the residual is obtained from equa- 
tions ( 4 ) and ( 5 ) of reference 1 or equation (26^) of reference 4. (in 
that eqmtion, the last term in the second parenthesis, 2cup sin a, 
should have a minus sign.) 


/ sin a ^ In T* \ 1 


T* W, 




'cp 


acp 


2U* 


sin a 
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TABIE I - CCMOTATIOII CF FLOW OH MEM sraEAMLIHE MD FIRST MD SBGCaiD 9-IKRIVM1IIVE3 AT MEAN BTRZAMLIHE - Continued 
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TABLE lU - EESrOUALS OBTAIKEID IN CEBTRIFUQAL IMPELLER 
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Station 1 2 3 4 5 6 7 8 9 10 12 14 16 18 20 23 

(t) Centrifugal ocanpresBor 


Plgure 2. - Concluded. Blade sections on mean surface of revolution. 
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(c) Eesultant. 

Figure 6. - Concluded. Comparison of velocities for present and relaxation 

solutions for turtine cascade. 
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(a) (T*r*p*W|)j. for totIoub stations. 

Figure 11, - Variations on chosen streamline for progressive cycles of ccmputatlon for centrifugal lameller. 
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(b) (P*W|)c for TarlouB stationa. 

Figure 11. - Cancluded. Yarlatlons on. clioeen atreamllne for progreBslTe cycles of ccanputation for 

centrifugal in^eller. 



NACA TN 2702 


43 



Figure 12. - Theoretical and experimental velocity distrlhution 
near "blade surfaces for centrifugal Inipeller. 
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Figure 13. - Theoretical and. experimental static-pressure dlstrihution 
along hlade surfaces for centrifugal inpeller. 
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